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Abstract 

The variational equation for the mean square displacement of the electron in the polaron worldline 
approach to quenched QED can be cast into a form which closely resembles the classical Abraham- 
Lorentz equation but without the conceptual and practical diseases of the latter. The connection 
with delay equations describing held retardation effects is also established. As applications we solve 
this integro-differential equation numerically for various values of the coupling constant and cut-off 
and re-derive the variational approximation to the anomalous mass dimension of the electron found 
recently. 



1 Introduction 

There is a long history of attempts at classical models of the electron until Quantum Electrodynamics 
(QED) took over as the most successful and precise theory in the microscopic domain [1]. Noteworthy 
among these models is, in particular, the one due to Abraham and Lorentz nearly one hundred years 
ago which is covered in many textbooks (see for example, Ref. [2]). It describes a classical electron 
under the influence of both an external force as well as the back reaction from the energy loss due to 
radiation of photons, and leads to a third-order differential equation which, however, has extraneous 
unphysical run-away solutions. The usual method to avoid these solutions is to convert the original 
equation into a second-order integro-differential equation. The classic nonrelativistic result is of the 
well-known form [2] 

^x(t) = -/ dse- s F(t + sr) , (1) 

dt A m Jo 

where r is the characteristic time 2e 2 /3mc 3 , m is the electron mass and F is an external force which 
can depend on time either explicitly or implicitly through x(t). The solutions to this equation do 
not exhibit run-away behaviour, however Eq. (1) clearly violates (nonrelativistic) causality because 
the motion at time t depends on forces at time t' > t. It is noteworthy that the analogous equation, 
without runaway solutions, for a relativistic Dirac particle [3] in 3 dimensions has only been derived 
rather recently [4]. 

From a modern perspective, of course, it is appropriate to wonder how nonrelativistic and classical 
equations like the one above are connected to a quantum mechanical and ultimately a quantum 
field theoretical description. In nonrelativistic QED this correspondence limit has been examined in 
particular by Moniz and Sharp [5] while Johnson and Hu have derived an Abraham-Lorentz-Dirac 
equation as a semi-classical limit within scalar quantum field theory [6]. 

In the present paper we discuss the fully relativistic field theoretic problem of a spin-1/2 electron 
dressed by an arbitrary number of photons in the quenched approximation to QED. Our approach 
is patterned after Feynman's celebrated variational treatment of the polaron problem [7], which was 
first applied by Mano [8] to a relativistic scalar field theory and re-discovered and expanded by us in a 
series of papers [9, 10]. Its main features are the description of relativistic particles by worldlines [11] 
parametrized by the proper time, an exact functional integration over the light fields (i.e. photons) 
and a variational approximation of the resulting effective action by a retarded quadratic trial action. 
In recent work we have extended this approach to more realistic theories, in particular to quenched 
QED [12, 13]. 

The approach is fully nonperturbative, Lorentz covariant, respects gauge symmetry and contains 
the exact one-loop self-energy in the small coupling limit. It does not rely on a perturbative, semi- 
classical or a derivative expansion. In Ref. [13] we concentrated on the divergence structure and 
renormalization of the theory, resulting in a compact expression for the anomalous mass dimension 
of the electron. In the present work we go beyond this and calculate also the finite contributions. 
We do this by numerically solving the variational equation for the position of the pole of the electron 
propagator for a range of coupling constants and UV cut-offs. We shall also show that the relevant 
variational equation can be written as an Abraham-Lorentz-like equation for the mean square displace- 
ment of the electron. As this "Variational Abraham-Lorentz Equation" (VALE) turns out to be (at 
most) a second-order integro-differential equation, no run-away solutions occur. Indeed, anticipating 



the result derived in the Section 3, in general this equation may be written in the suggestive form 

^(x\t)) = 2 (i 2 (t)) + ^-£™dsf(s){([Ax(t)} 2 ) - ([Ax(t - s)] 2 )} . (2) 

Here £„(£) is the Euclidean position at proper time t (t > 0), Ax M (i) is its deviation from a classical 
straight-line path, ke is a reparameterization parameter playing the role of a mass and < . . . > refers 
to a certain averaging over all worldlines occurring in the functional integral for the theory. These 
definitions will be made more precise in the following two Sections. At this stage we merely point out 
the marked similarity of this equation to Eq. (1), in particular if one re- writes the latter as 

cfi 9 r°° 

^x 2 (t) = 2± 2 (t) + — x(t)- / dse- s F(t + sr) . (3) 

dt z m Jo 

Of course, Eqs. (2) and (3) are not completely identical, nor should they be. In the absence of 
external forces the latter equation yields a differential equation with the solution that the particle 
moves uniformly along a straight line. On the other hand, there are no external forces in Eq. (2). The 
first term on the r.h.s. is the classical term while the integral characterizes an 'internal force' which 
has its origin in the constant emission and re-absorption of virtual photons. This term is sensitive 
only to the (random-walk like) deviations from a straight line. 

In the next Section we will describe the derivation of the variational equations, leading to the 
VALE (2) in Section 3 and will solve it numerically in Section 4. Understanding these numerical 
results will also give a new, simpler, derivation of the variational approximation to the anomalous 
mass dimension than the one presented in Ref. [13]. Approximate forms of the VALE are discussed 
in an Appendix. 

2 The variational method 

In order to introduce our notation and terminology, we begin with a brief summary of the essential 
points of the variational technique applied to quenched QED and refer the reader to Refs. [9, 12] and 
in particular Ref. [13] for the details. 

2.1 Worldline formulation 

The starting point is the worldline, rather than field theoretic, description for the propagator of 
a spinning particle which involves functional integrals over a bosonic worldline x«(£) as well as a 

Grassmannian Cu(t) which characterizes the spin. In Minkowski space time (with metric (+, )) 

the free part (modulo boundary terms) reads 



S = I dt 
Jo 



— x +iC,-C, + —x-C,x 



(4) 



where kq is the reparametrization parameter and x the supersymmetric (SUSY) counterpart to the 
proper time T. The dependence of the action on these two degrees of freedom is connected by a 
supersymmetric transformation [12]. Due to time-translational invariance only the time interval T 
matters and therefore in reality the integration limits in Eq. (4) and all following expressions are 



[to,T + to] with to arbitrary (e.g. or — T/2). The Gaussian functional integral over the photon field 
in the interaction part can be carried out analytically, resulting in a bi-local effective action 

Si = -^J^dhdt 2 |^A C iV (k)J„(k,x 1 ,Ci)M-k,x 2 , C2 ) e -i*-(si-* 2 ) (5) 

where G^ u (k) is the photon propagator in an arbitrary (covariant) gauge, 

2 
Jn(k,x,Q = &H C^fc-C (6) 

Kq 

is the (convection and spin) current of the electron and x\ = x(t\) etc. There is an elegant, compact 
and manifestly supersymmetric formulation in terms of "superpositions" and "superderivatives" [12] 
but in the following we will use the more cumbersome but explicit decomposition into bosonic and 
fermionic worldlines. 

2.2 Feynman-Jensen stationarity and variational principle 

Having integrated out the photons, the functional integrals over x and £ cannot be performed exactly 
but may be approximated via the introduction of a trial action St[x, (] and by making use of Feynman's 
variational technique: 

e iS = J e *s t .J_^_ = J e iS t . $ el{S ~^f St = Je iS < ■ (e l(S - St) ) St = Je iS * ■ e^ 3 ^ . (7) 

Even with Euclidean times, where exp(iS') — ► exp(— Se) , the usual Jensen inequality does not hold 
anymore since we are also integrating over Grassmann- valued trajectories. However, only stationarity 
of the above expression under variations of the trial action is required in the following. The subscript 
St on the average reminds us that the weight function is e lSt ; the averages in Equation (2) are to 
be understood in the same way. With an arbitrary trial action the above approximation becomes 
exact at the stationary point. In practise, however, only trial actions at most quadratic in x and £ 
can be used. We work with a trial action of this form obtained by making the modified free action 
Sq = Sq + p • x bi-local 1 



f T 
St = Ai p ■ x + / dt\ dt 2 

Jo 



T 

(8) 



--^9b{<j)xi ■ x 2 + i g' F {a) Ci ■ (2 - ij;crg' M (cr)xi ■ ( 2 X 



In addition to a scalar variational parameter Ai this trial action contains three arbitrary (but even) 
variational 'retardation' functions of a = t\ — £2 multiplying the three quadratic combinations x x (i.e. 
'bosonic', or B), £C ('fermionic', F) and x(, ('mixed', M). Only the first two of these are relevant 
in the calculation of the pole position of the electron propagator, which is obtained for T — ► 00. An 
explicitly supersymmetric trial action would require only one retardation function 

9b(o-) = gF{<?) = 9m(o-) = g{<?) ■ (9) 

As the exact action contains SUSY-breaking boundary terms (which will be discussed in more detail 
elsewhere [14]) we shall not enforce Eq. (9). It is therefore advisable to take the most general ansatz 

1 The integration over the endpoint x(T) — x of the trajectory with weight exp(ip ■ x) (where p is the external 
momentum of the particle) is included in the average. Note that the present retardation functions have a different 
normalization than in Ref. [101. 



and let the variational principle choose the optimal solution within the given class of test functions 
2 . Note that for Ai = 1, 5b(<t) = 9f(&) = 9m(&) = 5(a) the trial action (8) reduces to the free action 
(4). Because averaging with the free action is equivalent to first-order perturbation theory this implies 
that the variational approach gives the correct one-loop self energy for small coupling. 

2.3 Mano's equation 

The averages involved in the last term of Eq. (7) may be separated into averages of So — St (which, 
together with j exp(iSt), may be combined into a quantity 0) and the average of S± (denoted by V). 
The roles of O and V are not unlike those of, respectively, kinetic and potential contributions in a 
standard quantum mechanical variational calculation. Near the pole the electron propagator takes 
the form 

and the on-shell limit of the argument of the exponential of the last term in Eq. (7) directly yields the 
relationship between the electron's bare mass Mq and its physical mass M. Explicitly, this equation, 
termed 'Mano's equation', becomes 

M 2 = M 2 (2X-X 2 )-2ln[A B }-n[A F } + V[fi 2 B ,fi 2 F ] J . (11) 

Here the 'profile functions' A B and Ap are Fourier transforms of the bosonic and fermionic variational 
retardation functions, respectively: 

/+oo roo 

dagi(a)e iEa = 2 / da g t (a) cos (Ea) , i = B, F . (12) 

-oo JO 

They are fixed, as is A = Ai/Ab(0), through the Feynman- Jensen variational principle which guar- 
antees that Mano's equation (11) is stationary w.r.t. their variation. Once fixed, all quantities of 
interest in the field theory (e.g. masses, form factors, scattering cross sections [9]) may be expressed 
in terms of these. A B and Ap only appear implicitly in the potential V, through the bosonic and 
fermionic 'pseudotimes' fj, B F (a) defined by 

M-tr^wm^ifY (13) 

2.4 Kinetic and potential terms 

Explicitly, the kinetic terms Q[Ai] are given, in 4-dimensional Euclidean space, by 

"M = ^fM'°^) + ^ - i) • (") 

Here k b > is a parameter which re-parameterizes the proper time without changing the physics. It 
is useful to keep it because it plays the role of a mass in the worldline description (see Footnote 3 in 
Ref. [10]; the Euclidean parameter up is related to its Minkowski counterpart kq through kq = ike)- 

2 It should be noted that Eq. (8) is still not the most general quadratic ansatz as the retardation functions are 
assumed to depend only on the proper time difference and additional Lorentz structures are absent. For the general case 
in a scalar theory see Ref. [10], Appendix C. 



The potential term V results from averaging the bi-local effective action (5) and performing the 
limit T — ► oo. While details of this calculations are rather involved and will be given elsewhere [14] 
the final result was already presented in Ref. [13]. As V is ultraviolet (UV) divergent we write it down 
in d = 4 — 2e dimensions. In Euclidean time (with p E = —M 2 ) it takes the form 



V 

V-i 

v 2 



V1 + V2 

v (d-1) / da - — -T — : ;: 

k e Jo J (2vr) d k 2 + m 2 

d d k 1 



^^X 2 Tda 
he Jo 



(2vr) d k 2 + m 2 



M 



Ji F {a) 
2 (d-2) 



Ab(^ 



■VE? 



k 2 



E(a,k) 
E(a,k) . 



Here v is an arbitrary mass parameter, the abbreviation 



E(k,a) 



exp 



1 



2k e 



k jU B (cr) — 2A/c • ps a 



(15) 
(16) 

(17) 



(18) 



has been used and a photon mass m has been kept in the photon propagator. The separation of V 
into two parts makes sense because of several facts: first, it is seen that V2 is more singular than V±, 
having an additional power k 2 in the integrand. Second, V2 vanishes for massless electrons M = 0, is 
repulsive (V2 > 0) and only depends on the bosonic pseudotime. 

Previously the anomalous mass dimension jm was calculated analytically in dimensional regular- 
ization [13]. However, from a numerical point of view this regularization is extremely cumbersome in 
nonperturbative calculations (see, for example, the Dyson- Schwinger equation studies in Ref. [15]) and 
is usually replaced by a simple momentum cut-off. Proper time regularization, i.e. a lower cut-off of 
proper time integrals at a = 1/A 2 , A —> 00 could also be used (it would maintain translational invari- 
ance), but it violates reparametrization invariance. In the present work, it is much more convenient 
to make use of a form factor which is very similar to the nonlocal regularization method proposed 
in Ref. [16]. Here the main reason is that the UV divergence in V originates in an undamped (for 
a —> 0, in which case jUj((j) — ► a) Euclidean momentum integral over the exponential factor E(k,a). 
This is the only appearance of [X 2 B in V and hence, after introducing a form factor of the form 



F(k 2 



exp 



2A 2 



(19) 



all dimensionally regulated expressions in Ref. [13] may be converted to form factor regulated ones 
by setting d = 4 and making the simple replacement 



^K CT ) 



k e 



As(°") = Mb(^) + ^2 



(20) 



With this regularization understood, the different terms in the potential V of Eq. (15) become (after 
performing the momentum integration for m = 0) 



ViIpbiPfI 



v^l] 



3a < x r/.2f~Vk2 (.:2 f „\\2 



47T Jo 



AbIo-J 



3a 



7T 



-KE 



da 1 - [1 +7(cr)]e' 



-7(0-) 



cr 



7(a) 



(21) 



(22) 



Here 



7(a) 



X 2 M 2 a 2 



fi 2 (a) 



dn 2 (a) 



(23) 



2ke fi%(cr) ' da 

and a = e 2 /(47r) (in the real world ~ 1/137 ) is the fine structure constant. 

It is clear from Eqs. (21,22) that the a — ► UV-divergences due to inverse powers of the 
bosonic pseudotime are now regulated and that both V\ and V2 as well as fij do not depend on 
the reparametrization parameter ke, which only sets the scale for the a and E- variables 3 . One 
also sees that bosonic and fermionic degrees of freedom do not enter symmetrically in the interaction 
(e.g. V2 only depends on jj? B ) and therefore the pseudotimes are in general different. This is due to 
supersymmetry violation by boundary terms in the trial action and could affect the spin structure 
of the propagator in our variational scheme. Since the spin-dependent terms are subleading in the 
limit T — > 00 and are not considered in the present investigation we just let the variational principle 
decide how much of SUSY violations it tolerates with the present trial action. Fortunately, as we will 
observe numerically and is to be expected from boundary terms, these supersymmetry violations are 
restricted to large values of a (or small values of E) and therefore do not have any influence on the 
UV behaviour of the solutions. 

Because the fermionic contributions, both in the kinetic term of Mano's equation (11) as well as 
in V\ (Eq. (21)), appear with an opposite sign to the bosonic contributions they would cancel them 
for exact supersymmetry, i.e. fip(a) = [i 2 B (a). Note that for a — > a cancellation of this sort is 
absolutely necessary as otherwise one would encounter the same quadratic divergences which occur 
in the propagator of scalar QED. In the context of the variational calculation the cancellation occurs 
only after variation because one needs a restoring force in the variational principle. A simple example 
illustrating this subtlety is given in Ref. [17]. 



2.5 Variational equations 

With the explicit expressions for 0, and V above it is now a simple, albeit somewhat tedious, exercise 
to derive the variational equations for A, Ap and As- One obtains, respectively, 



A 



3a ke 

2^XM 2 



da 



(^) 2 



(A 2 ^ 2 



fi% 



B^~ e -, + 4(l + 7 + 7 2 )e^-l 



1 



_2_ f™ j 5V sm 2 {Ea/2) 

~ E h M JIW) E 2 A t (E) 



(24) 



(25) 



ME) 

where we define (— ) = 1 and (— ) = —1 , the difference in sign originating from the opposite sign 
of OLAb] and fif-Ap] in Mano's equation. The functional derivatives in Eq. (25) are 



5V 



3o 



d 



5fi 2 F (a) 



— k e ~r 
2ir da 



£Hg) e -7W 



and 



5V 



Sn 2 B {a) 



3a d 

2ir da 



/4(°") 



3a 1 
Ms 



4tt 



(/i 



2 \2 
F> 



(A? 



Ab 



- (2 - 7)e"T 

4 (l + 7 + 7 2 )e~^-r 

j2 Zj 



^This is due to the reparametrization dependence of Ai(E, ke) — Ai(keE) and jj% (cr, ke) — k-e ^{ct/ke) 



(26) 



(27) 



In the limit in which the interactions are turned off, we have Ai(E) = A = 1 and therefore fJ%(a) = a, 
fij = 0, M = Mq. In the interacting theory the pseudotimes fJ^(a) still have linear behaviour for both 
a — > 0, oo. From Eqs. (26,27) we then see that without regularization 



SV\ const. SVo const. 



6 ft? (a) a 3 Sfi 2 B (a) a 



(28) 



for a — ► 0. It is this more singular UV-behaviour of SV\ which requires regularization also in the 
variational equations. This is the trademark of a renormalizable theory whereas the more benign 
behaviour of 6V2 is characteristic for a super-renormalizable theory. On the other hand in the infrared 
region V2 dominates because the very last term in Eq. (27) is not exponentially suppressed 

5V 2 ct-^oo const 
6fi%(a) a 4 

This is due to taking a vanishing photon mass m = in Eq. (17). 

3 The variational Abraham-Lorentz like equations 

It is obvious that the profile functions Ai and the pseudotimes fif contain the same information, as 
seen in Eq. (13). It would be more efficient, therefore, to eliminate one of these from the variational 
equations (25). We shall now show that it is possible to rewrite Eq. (25) entirely in terms of the 
pseudotimes. 

A differentiation of Eq. (13) with respect to a yields 

-i(s d 1 1 , 2 f°° jtp s[n ( E(7 ) ( ! 



^-■S^ + Hk iE ^\ME)-')- (30) 

It is important here to realize that the pseudotime is even (as can be seen from Eq. (13)) and that 
the first term on the r.h.s. therefore gives sgn (a). The remaining integral is convergent due to 
subtraction of the asymptotic value l/Aj(oo) = 1. Another differentiation leads to 

2 f 00 / 1 \ 

/i?((7) = 2 6(a) + -J o dE cos(Ea) {j— - lj . (31) 

If we now insert Eq. (25) we obtain 

2 1 r°° r°° ■ 5V 1 r i 

fi(a) = 26(a) / dE da > (-y-^— 1 - cos(Ea') cos(Ea) . (32) 

vr k e Jo Jo <Vf (°" ) E z Ai(E) L J 



By using the addition theorem for the cosine function and the definition (13) the ^-integration can 
be performed exactly and gives 



1 f°° , , , w SV 



1 f 00 
fi (a)-— da'(-) 

KE JO 



M?(<r) - ^ ( a + a ') " 2^ ^ ~ a '' 



2 6(a) , a >0 . (33) 



ke Jo S^H* 7 ' 

The appropriate boundary conditions for solutions to this integro-differential equation are 

M ?(0) = 0, lim J4(a) = 1. (34) 



The proper times a, a' are restricted to be nonnegative and we therefore have to take the absolute 
value in the argument of the shifted pseudotime in Eq. (33). Alternatively, this restriction may be 
avoided by remembering that /j, 2 (a) is an even function of a and hence the integrands appearing in 
the expression for V (Eqs. (21,22)) are also. One obtains 

1 f +oc ■ 5V r i 

£*(*) - 7T- / da' (-Y -^- tf{a) - tf{o - a') = 2 5(a) , -oo < a < +oo (35) 



2ke J-oo fy^V) 

and, in fact, even the 5-function may be eliminated by performing differentiation with respect to |<r| 
rather than a: 

1 f+°° f)V r i 

^"(M)-*— / da* (-y-^— U(\a\)-tf(\a-a'\)] = 0, -oo < a < +oo (36) 

2k E J-oo d/j,f(\a'\) L J 

Here " denotes differentiation with respect to the argument. 

We shall present numerical solutions to these variational equations (which, as promised, no longer 
involve the profile functions directly) in the next Section. At this stage, however, we would first like 
to discuss the meaning of Eq. (36) for the bosonic case, i.e. i = B. In Ref. [10], Eqs. (25,26) the 
expectation value of x^(ti) — x^fa) and [x M (£i) — x^fo)] 2 , when averaged with the trial action, were 
calculated. The corresponding results for the present case (i.e. QED, with Euclidean metric, and 
setting £2 = and hence x^fo) = for convenience) become 

(HVDst = f = kdasB.(*)l 

he 

(x 2 (t)\ = (\x(t)\)] l+±^(t) = xL M .(t) + ([Ax(t)] 2 \ , (37) 

i.e. the mean squared displacement of the electron is made up by an overall quadratic drift due to the 
electron's momentum, i.e. growing like t 2 , and a term proportional to the pseudotime, growing like 
t for small t; the latter characterizes both the quantum mechanical Brownian motion as well as the 
continual random "kicks" from emission and absorption of virtual photons in the cloud surrounding 
the bare particle. 

Inserting Eq. (37) into the variational equation (36) results in the Abraham- Lorentz-like equation 
mentioned in the Introduction, i.e. 

5V/5fj> 2 may then be interpreted as (internal) "force" acting on the bare electron. 

As mentioned in the Introduction, the original Abraham-Lorentz equation actually is a third-order 
linear differential equation for the position of the electron rather than an integro-differential equation. 
It is possible to write the VALEs in Eqs. (33) - (36) in a similar way by converting them into a special 
form of delay equations. This class of differential equations has been studied extensively for problems 
of radiation damping in electrodynamics and general relativity and in mathematics [18]. Indeed, by 
invoking the first mean value theorem (this requires that 5V/5fj, 2 (a) is of one sign, which appears to 
be fulfilled) to the integral in Eq. (36) one can replace the <r'-argument of the pseudotime by some 
mean time < r% < 00. The remaining integral then just gives a constant 

r°° ■ SV 

C\ = / dai-y-^- (39) 

Jo (>i2f(a) 



Therefore we have for a > 

1 . r _ 1 - i - i 

= 0. (40) 



fi{a) - —C\ 



UK*) - 2$ (°" + T i) - ^ (l ff ~ Ti ' 



K E 

The above (seemingly linear) delay equation is exactly equivalent to the original (nonlinear) VALE 
because the constant C\ is in general a functional of the pseudotime and because by construction the 
delay T{ will also depend on the external proper time a as well as on fsj(cr). It is worth noticing that 
C\ also governs the asymptotic behaviour of the profile function 

a ^°°> = 1+ oHf + "- < 41 > 

By differentiating Eq. (36) it is seen that it also determines the initial "jerk" [19]: 

t (a - +0) = ^ . (42) 

K E 

For the fermionic case (or in the supersymmetric limit) 5V/d/j, F is a total derivative (see Eq. (26)), 
so here one finds the exact expression Of = 3aA 4 /(2-7r). 

In many applications of delay equations a constant delay is assumed. Although this does not seem 
to be a valid approximation for QED it is clear that the delay n is a remnant of the photon degrees 
of freedom which have been integrated out and no longer appear in the equation of motion for the 
mean displacement of the bare electron. Since in QED the "force" 5V/5/J, 2 is very singular at small 
a one expects a very small delay r ~ 1/A 2 (all integrals are dominated by small values of a'): the 
average time a virtual photon is "in the air" is very short, but the kick given to the bare electron is 
very violent. 

Finally, we emphasize that in the discussion above no use has been made of the specific form of 
the interaction, i.e. the VALEs (33) - (36) hold for any dressed particle. If we concentrate now on the 
specific form (21,22) of the interaction in QED further simplifications can be made. For investigations 
of the divergence structure of the theory as in Ref. [13] it should be a reasonable approximation 
to neglect Vi since it it was shown to be less singular than V\. Then 5V/5/j, 2 (a) is just a total 
derivative (see Eq. (26)) and one may integrate the VALE with respect to a. In this supersymmetric 
approximation fermionic and bosonic pseudotimes are identical and for a > obey the equation 

3« f°° , / AJUSYJO -y(a') 



4vr Jo [^susyKJJ l j 

Since the regularized pseudotime just adds a constant to the pseudotime // 2 (<r), Eq. (43) also holds 
for /i 2 (o"). To exhibit the UV divergences one may even neglect the electron mass altogether. In this 
"asymptotic SUSY" (ASUSY) approximation the quantity 7 in Eq. (43) is set to zero which allows 
a further integration over a. Taking into account the boundary condition fJ- 2 (0) = one obtains the 
simpler integral equation (a > 0) 

„2 (^ , 3a f°° j> ^ASUSY (g + gO + ^ASUSY ( W ~ °' I ) ~ 2 MJsUSY (O _ „ (AA ^ 

MasusyIo-J + — / da r^ 7-77 - a. (44) 

4vr Jo ^asusyI^ 7 ) 

Adding ke/A 2 on both sides turns it into a compact, but still nonlinear integral equation for 

( ° k2 \ A ' ~2 

V\8 = : = — ^asusyW ( 45 ) 

\ KE J K E 

10 



namely, 

y(s + af) 



3a r + °° 
y(s) = 1 + \s\ - — / ds' 

27T J-oo 



*/(*') 



1 



(46) 



We have been unable to find an exact analytical solution for this scale- free equation. For the purpose of 
calculating the anomalous mass dimension it is, however, sufficient to know the asymptotic behaviour 
of y(s) which will be derived below. An approximate, third-order differential equation (i.e. of the 
same type as the original Abraham-Lorentz equation) is obtained in Appendix A. 

4 Numerical results 

Without any approximation the three variational equations, i.e. the fermionic and bosonic VALEs 
(Eq. (33)) together with the variational equation for A (Eq. (24)) may be solved iteratively by stepwise 
numerical integration over a. Because the profile functions Ap^ have been eliminated in favour of 
the pseudotimes, it is not necessary to numerically evaluate the cosine transform (Eq. (13)) for the 
pseudotime and its derivative. This simplifies the numerical effort as compared to the published 
variational treatments of super-renormalizable and finite theories [9, 10, 20]. 

However, the presence of two vastly different scales characterising both the ultraviolet ( a ~ 
ke/^ 2 ) and the infrared region ( a ~ ke/{XM 2 ) with A — ► for A — * oo, see below) makes the 
calculation much harder for the present renormalizable theory. Stable numerical results have been 
obtained for cut-offs as large as A = 500 M and coupling constants up to a = 0.5 — 1.0 and have been 
tabulated in Table 1. The cut-off needed to be decreased somewhat as the couplings were increased 
in order to retain numerical stability. As in previous work (cf. Refs. [10, 20]), we have checked the 
accuracy of our results by comparing the direct evaluation of the kinetic terms f2j with those obtained 
with the help of a virial theorem which relates the ilj to the potential V through the use of the 
variational equations: 



a 



(-)' / da [i4(<?)-^M)\ J^-y (47) 



A relative accuracy of order 10 -5 was achieved for not too large cut-offs. 

The results in the Table clearly demonstrate, as anticipated in Sec. 2.4, that supersymmetry is 
almost perfect. This may also be seen in Fig. 1 where the pseudotimes have been plotted for one 
value of a and A. Differences between the bosonic and fermionic pseudotimes only show up in the 
infrared region where the less-singular potential contribution V2 dominates. This also leads to slightly 
different values of the bosonic and fermionic profile function at E = 0. Our previous conjecture that 
SUSY violations would vanish with A — > 00 [17] turns out to be unfounded: in Table 1 we also give 
the quantity 

which may be considered as a measure of the importance of SUSY violations in Mq . One sees that 
for the whole range of coupling constants which we consider they remain at the percent-level or below 
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a 


k/M 


A 


A B (0) 


A F (0) 


M /M 


10 3 A 5 


M 2 \ 


/3(a) 


XA B (0) 


0.1 


50 


0.69862 


1.4009 


1.3992 


0.81663 


-0.36 


0.95457 


0.95450 


0.97869 




100 


0.65402 


1.4964 


1.4946 


0.79011 


-0.37 


0.95452 




0.97866 




300 


0.58901 


1.6615 


1.6595 


0.74981 


-0.37 


0.95451 




0.97865 




500 


0.56102 


1.7444 


1.7423 


0.73178 


-0.37 


0.95450 




0.97865 




1000* 


0.52515 


1.8636 


1.8614 


0.70800 


-0.34 


0.9545 




0.97865 


0.2 


50 


0.48548 


1.9802 


1.9716 


0.66592 


-1.31 


0.91342 


0.91332 


0.96132 




100 


0.42570 


2.2582 


2.2484 


0.62355 


-1.38 


0.91335 




0.96129 




300 


0.34559 


2.7815 


2.7696 


0.56182 


-1.38 


0.91333 




0.96128 




500* 


0.31366 


3.0647 


3.0516 


0.53522 


-1.42 


0.9133 




0.96128 


0.3 


50 


0.33610 


2.8171 


2.7923 


0.54267 


-2.85 


0.87619 


0.87610 


0.94682 




100 


0.27631 


3.4266 


3.3966 


0.49201 


-2.88 


0.87612 




0.94681 




300 


0.20252 


4.6752 


4.6342 


0.42122 


-2.88 


0.87610 




0.94680 




500* 


0.17527 


5.402 


5.355 


0.39184 


-2.99 


0.8760 




0.94680 


0.4 


50 


0.23214 


4.0256 


3.9682 


0.44224 


-4.68 


0.84249 


0.84241 


0.93449 




100 


0.17914 


5.2165 


5.1423 


0.38847 


-4.71 


0.84242 




0.93448 




300* 


0.11877 


7.868 


7.756 


0.31629 


-4.78 


0.8423 




0.93448 




500* 


0.09810 


9.525 


9.390 


0.28746 


-4.81 


0.8423 




0.93448 


0.5 


50 


0.16015 


5.769 


5.651 


0.36059 


-6.75 


0.81191 


0.81186 


0.92384 




100 


0.11616 


7.953 


7.791 


0.30709 


-6.77 


0.81187 




0.92383 




300* 


0.06981 


13.23 


12.96 


0.23805 


-6.88 


0.8118 




0.92383 


0.6 


50 


0.11046 


8.828 


8.056 


0.29430 


-8.95 


0.78410 


0.78406 


0.91451 




100 


0.07541 


12.127 


11.800 


0.24315 


-8.98 


0.78406 




0.91450 




300* 


0.04117 


22.21 


21.61 


0.17966 


-9.06 


0.7840 




0.91450 


0.7 


50* 


0.07623 


11.890 


11.485 


0.24048 


-11.2 


0.75867 


0.75869 


0.90626 




100* 


0.04905 


18.48 


17.85 


0.19290 


-11.2 


0.75867 




0.90625 


0.8 


50* 


0.05266 


17.07 


16.37 


0.19680 


-13.6 


0.73540 


0.73544 


0.89889 




100* 


0.03198 


28.10 


26.95 


0.15336 


-13.6 


0.73535 




0.89887 


0.9 


50* 


0.03644 


24.49 


23.30 


0.16130 


-15.9 


0.7140 


0.71407 


0.89227 




100* 


0.02091 


42.66 


40.60 


0.12219 


-15.9 


0.7139 




0.89225 


1.0 


50* 


0.02526 


35.08 


33.13 


0.13244 


-18.2 


0.6943 


0.69435 


0.88627 




100* 


0.01372 


64.61 


61.01 


0.09757 


-18.3 


0.6941 




0.88622 



Table 1: Some results of the numerical solution of the variational equations with form factor regu- 
larization for different cut-offs and coupling constants. A relative accuracy of 10 -6 was required at 
each point between subsequent iterations; in entries marked by a * only 10™ 5 relative accuracy was 
achieved. The quantity As measures the SUSY violation and is defined in Eq. (48). The last three 
columns compare the numerical results with analytical approximations explained in the text. 
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Figure 1: The pseudotime divided by a (solid line) as a function of a (on a logarithmic scale) for 
coupling constant a = 0.7 and cut-off A = 100 M. Bosonic and fermionic pseudotime are indistin- 
guishable in this plot. The dashed curve shows the difference of both. 

and become cut-off independent within numerical accuracy. 

The ratio a ft 2 (a) / fj 2 (a) is plotted in Fig. 2 for a = 0.5 and various cut-offs. Independently of 
the actual magnitude of /U 2 (u), power-like behaviour of the pseudotime corresponds to a horizontal 
line on this plot. Because /i 2 (c) — ► a and ^ 2 (<r) — ► a/A(0) all curves go towards unity for small 
and large a. What is interesting, however, is that for an increasing range of intermediate a-values (as 
A — > oo) the pseudotimes exhibit power-like behaviour of the form cr@, [3 < 1. Moreover, the value of 
(3 seems to become independent of the cut-off as A increases. 

This result can be understood analytically: for large cut-offs and moderate values of a the ap- 
proximate VALEs (43) and (44) should become valid. Inserting the ansatz 



/i 2 (a) 



s -a 



ft 



into Eq. (44) one obtains 
s <7 /3 + 



3a 

47T 



da 



1 + 



a" 



+ 



1 



a 



a 1 



a + 



ke_ 
A 2 ' 



(49) 



(50) 



Scaling a' = to shows that the integral is proportional to a and therefore dominates the l.h.s. for 
large values of a provided f3 < 1 which turns out to be the case for positive couplings . Hence, for 
large a, the power (3 is only a function of the coupling constant a and determined by the implicit 
equation 

t 3a 



dl 



lx/3 

1 + 7» + 



(51) 



4 For a < no power- like solutions exist anymore as the term o^ would dominate for /3 > 1. 
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0.95 



0.90 



0.85 
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Figure 2: The effective power of a for the the bosonic pseudotime n 2 (cr) from numerical solutions of 
the variational Abraham-Lorentz equation for QED. As in Fig. 1, the difference between lSp(cr) and 
/ u s( <T ) ^ s too small to show up on this plot. Results for various cut-offs A in units of the electron mass 
M are shown. The dashed horizontal line denotes the numerical value from the analytical solution 
(52) at coupling constant a = 0.5. 



After performing the integral one obtains the following transcendental equation for (3(a) 

ft r, [ft r\ 27T 

2^ tan 2^ = Ta ■ 



(52) 



Numerical solutions for a variety of a's have been tabulated in Table 1. For small coupling constants 
the solution to this equation behaves like 



3a [3a\ 

"^-aF + UF.I + 



(53) 



confirming that for positive coupling constant (3 < 1, while for large a it goes like (3 — ► ^/8/(37rct) — v 0. 
Fig. 2 shows that the value (3(0.5) = 0.81186 (see Table 1) explains the intermediate-range behaviour 
of the pseudotime very well. For very large a the ASUSY approximation (44) becomes invalid and 
the pseudotime reverts to linear behaviour. 



5 The anomalous mass dimension 

The power- like behaviour of /x 2 (<r) is also the key for obtaining the variational approximation for the 
anomalous mass dimension regularized via a form factor rather than dimensionally [13]. To do this 
we first note that in V the variational parameter A occurs exclusively in the combination AM. Let 
us, therefore, define the alternative dimensionless variational parameter x\ = AM/A so that Mano's 
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equation becomes 

_ r m 

(54) 



Ml = A 2 



2x A — - x\ - 2H(a,x x ) 



where H(a,x\) is the dimensionless combination (Qb — &F + V)/A 2 . Note that, by construction, 

H(a,x\) no longer depends on M explicitly. Therefore, as it is dimensionless, it can also no longer 

depend on A explicitly. The dependence on the cut-off can only enter implicitly through the variational 

parameters and functions. However, the r.h.s. of Eq. (54) is stationary with respect to variational 

parameters/functions, so that the dependence of Mq on A, at fixed M, is just given by the explicit A 

dependence in Eq. (54), i.e. 

M 2 
2— 2--2Mza. (55) 



9M 2 



M fixed 

Writing Mq = Zm M v (where v is an arbitrary mass parameter) , the anomalous mass dimension may 
now be evaluated with the help of this flow equation as 



dlogZ M _ dlogMp _ _A_ Ml_ M \ _ M_ 2 
dlogi/ " Slog A ~ 2M 2 I A Xx J M 2 



1M= ^^ = -^^ = -^772 \ 2 ^r- 2M ^ =172 A " 1 - (56) 



Of course, the (existence of the) limit A — ► oo is understood in the above equations. If this limit 

exists then 7m must necessarily be cut-off independent so that a further differentiation of Eq. (56) 

with respect to A gives 

M 2 dM 2 M 2 <9A 

" ~m*^a x + mI M- (57) 

Using Eqs. (55) and (56) we obtain hd\/dh = —2jm A and integration then shows how the param- 
eter A behaves for very large cut-offs 

A A -Z!°? const . (^y M . (58) 

This behaviour comes as no surprise since we know [10] that in the worldline formalism the bare and 
effective mass of the quantum mechanical particle are ke and k^/A, respectively, for which a similar 
relation as for the actual masses is expected. Thus the anomalous mass dimension can be determined 
either numerically by solving the variational equations and evaluating Eq. (56) for larger and larger 
cut-offs or analytically from Eq. (58) by finding the cut-off dependence of A. 

In the following we will pursue the latter option which is possible as we know the approximate 
behaviour of the pseudotime for small, intermediate and large values of a (here we take for simplicity 

K E = 1) 



/2 2 (ct) ~ (a + ^ J 9 (<7i - a) + soa? @ {a - a ± ) @ (a 2 - a) + -£- 9 {a - a 2 ) 



(59) 



The regions are separated by o\ = xi/A 2 , where the cut-off 1/A 2 becomes effective in /2 2 , and 
o"2 = 2x2/(A 2 M 2 ^4(0)) where the exponential exp(— 7) becomes important and the ASUSY approx- 
imation breaks down. X\,X2 are numbers of order one which can be determined (together with the 
normalization factors ^4(0) and sq) by matching the approximate solutions at the boundaries. As we 
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only need the leading terms in the cut-off A in order to determine 7m, the actual numerical values of 
xi,X2 do not matter. Matching gives sq ~ (A 2 )' 3-1 as expected from dimensional arguments and 



_i 



/ A 2 \H 

A(0) ~ (j^j ■ (60) 

We now insert these expressions into the variational equation for A in SUSY approximation 

1 3a f°° j 1 (l+7 + 7 2 )e"^-l 

A 7r Jo A 2 ^) 7 2 



- = 1 + -L d °im^ J • ( 61 ) 



As 7 only becomes large in the region a > o"2 we may expand the exponential in the intervals for 
small and medium values of a. After performing the integrations, one obtains 



1 3a 



-ln{l + Xl ) 



1 
+ 



2(1 ~P) 



A(0) " 



1+Xl 



+ L4(0) ■ const] \ . (62) 



If A(0) would decrease with increasing A (or stay constant) this equation would imply that A would 
go to a constant. This, however, is inconsistent with Eq. (60). Hence A(0) must increase with A, 
with Eq. (62) implying that A^4(0) remains constant. This behaviour is clearly seen in the numerical 
results shown in the last column of Table 1. Also, in this case the low-a region can be neglected, 
which is reasonable because otherwise the final result would depend on the details of the regularization. 
Comparing this with Eqs. (58) and (60) one sees that consistency requires that 



1 " (M 

1M = —g— (63) 



and so Eq. (56) implies that 



MS 



M»A-"- (64) 

Comparison of the 8 and 9 columns in Table 1 shows clear numerical evidence that this relation is 
indeed fulfilled. Note also that Eqs. (52) and (63) imply that jm is a solution of the implicit equation 

-a = (1 + 7m) tan ( - M ) , (65) 

4 V 2 1 + 7m / 

in agreement with the result obtained in dimensional regularization [13]. The convergence of the 
perturbative expansion and the analytic properties of the solutions of such transcendental equations 
have been studied in Ref. [21]. The present derivation adds the insight that for a < no power-like 
solution of the VALE (44) exists and therefore Eq. (52) is not applicable anymore. Obviously, this 
peculiar behaviour when changing the sign of the coupling constant is not contained in Eq. (65) 
but in agreement with Dyson's old qualitative argument [22] that a = is an essential singularity in 
QED. We also note that Dyson- Schwinger calculations with a particular ansatz for the electron-photon 
vertex lead to a similar implicit equation [23] but a gauge dependence remains and solutions exist 
only below a critical coupling. 
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6 Conclusions and outlook 

In summary, we have shown that the variational formulation of worldline QED very naturally leads 
to an equation which is similar to the one considered much earlier by Abraham, Lorentz and Dirac in 
attempts to describe the electron and its self-interaction with the radiation field. In contrast to these 
attempts our approach contains (almost) all the ingredients of the relativistic field theory of electrons 
and photons, in particular its divergence structure. This has been demonstrated by numerically 
solving the variational Abraham-Lorentz equation (VALE) for a variety of cut-offs and by deriving 
an approximate nonperturbative expression for the anomalous mass dimension of the electron. We 
have shown how the approach leads naturally to qualitatively different behaviour of the theory for 
a > and a < 0. Furthermore, while the present investigation has been restricted to a free electron 
interacting with its own radiation field, the extension to the case when an external field is present 
as well should be straight-forward. This would allow a study of how this field-theoretically based 
worldline variational approach to QED avoids the pitfalls of pre-acceleration and acausality which 
have plagued all classical attempts. 

However, for further progress we deem it more important to first study the issue of supersymmetry 
(breaking) within the worldline variational approach in more detail. For massive electrons small 
violations of worldline supersymmetry have been observed which are "soft" in the sense that the 
ultraviolet behaviour of the theory is not affected. These violations presumably reflect the SUSY- 
breaking generated by the different boundary conditions for bosonic and fermionic variables in the 
exact action. Further investigation is required into what role, if any, this SUSY-breaking plays in the 
full spin structure of the electron propagator and how it is manifested in the worldline variational 
approximation. Such an understanding is required for future applications of this non-perturbative 
approach to physical processes. Finally it is not inconceivable that similar VALE's as derived here in 
the variational approach for the propagator will also emerge for the full interacting vertex. 



Appendix 



A Approximate 3 -order differential VALE 

Here we show how a 3 -order differential equation for the mean square displacement arises approxi- 
mately in QED. For this purpose we invert Eq. (13) to obtain 

1 l'°° 

ME) =1 + J du ti( u ) C0 < Eu ) ( A -!) 
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and observe that for small coupling the profile function stays close to unity. Therefore one has in this 
case 



a 



dE 



o 



2k e 

IT 

^TdE 
vr Jo 



1 



log ( 1 + 
1 



1 



ME) 



ME) 

2 

+ 0(1 



1 + 



1 



ME) 

1 ^ 3 



1 



ME) 



2 Jo 



#(*) 



Variation first gives a fourth-order ordinary differential equation 



d A 



sv 



do*™"' ' 6tf(a) 



0. 



(A.2) 



(A.3) 



but in the supersymmetric limit of QED the functional derivative of the "potential" V is a total 
derivative (see Eq. (26)) so that an integration yields 

■\2/ 



/./ 



3a Mo) , (ff) 
V 7 2vr /x 4 (cr) 



(A-4) 



This 3 rd order differential equation clearly is specific for QED and therefore it is tempting to interpret 
the "kinetic" term Q in this case as total (integrated) "radiative energy loss". Indeed, it is positive 
and in the approximate form (A.2) the integrand has precisely the form of the Larmor power formula 
(see, e.g. Eq. (17.6) in Ref. [2]). Without external force there should be, of course, no real radiation 
loss - in fact, for exact supersymmetry Qb is completely cancelled by ftp- Of course, the 4 order 
equation also applies (approximately) to other theories where this interpretation does not make sense. 

Since we have assumed /i 2 as small in Eq. (A.l) the present approximation should be similar to a 
derivative expansion. Indeed, if in the integrand of the VALE (43) 



fi 2 (a + a') - i? (\a - a'\) = fi 2 (cr> + a<) - fi 2 (cr> - a<) ~ /i 2 (o>) a < 
is expanded to first order and the resulting equation divided by a one obtains 



A 2 (a) 



3 a 
27 Jo 



1 + ^ / da' a' 



r,4 - ' ' 



/*V) 



+ 



3 a 

2^ 



'da'^e^ 



J&V) 



(A.5) 



(A.6) 



This can be converted to a differential equation by differentiating w.r.t. a and replacing the square 
bracket by l/(/i 2 — ajl 2 ). Another differentiation then gives 



•••2 . . 3a 
, (a) * - 



A 2 (a) 



a ft 2 (a) 



~^(a) 



(A.7) 



which is identical with the third-order equation (A. 4) if the factor in square brackets is replaced by 
its perturbative value 1. The boundary conditions are /U 2 (0) = 0, /t 2 (0) = 1 and a peculiar value of 
the second derivative at the origin 

2 



M 2 (0) = — / dafi 2 F {a 

ke Jo "ff 



6h%(<t) 



SUSY 



3a 

27 Jo 



da 



TO 



3 -7(ff) 



(A.8) 



expressed as functional of the still unknown solution. This follows from Eq. (36) in the limit a — > 0, 
together with Eq. (26) and can be considered as a reminder that the nonlocal actions we are working 
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with originate from local actions with the photonic degrees of freedom integrated out. In this way 
Ostrogradski's no-go theorem for higher derivative and nonlocal theories [24] is evaded. 

It is interesting that the (approximate) equation (A. 7) also allows for power-like solutions in the 
massless case whereas the perturbative version (A. 4) does not. Indeed, inserting the ansatz (49) into 
Eq. (A. 7) shows that the powers on both sides of the equation match and leads to the cubic equation 

(3-1 ~ ^(3 2 ((3-2). (A.9) 

This gives the correct limit (53) for a — > and the same l/y^-behaviour for large a as the exact 
result (52) only with a coefficient ir/^/8 = 1.11 times larger. That the derivative expansion works 
rather well also at intermediate coupling constants is demonstrated, for example, by the numerical 
value (3(0.5) ~ 0.81277 obtained from Eq. (A.9) to be compared with the correct (3(0.5) = 0.81186. 



References 

[1] S. S. Schweber: QED and the Men who made it, Princeton University Press (1994) . 

[2] J. D. Jackson: Classical Electrodynamics, 2nd edition, John Wiley (1975), chapter 17 

[3] P. A. M. Dirac, Proc. Roy. Soc. A 167 (1938) 148 . 

[4] M. Ibison and H. E. Puthoff, J. Phys. A 34 (2001) 3421 . 

[5] H. Levine, E. J. Moniz and D. H. Sharp, Am. J. Phys. 45 (1977) 75 ; 
E. J. Moniz and D. H. Sharp, Phys. Rev. D 15 (1977) 2850 . 

[6] P. R. Johnson and B. L. Hu, Phys. Rev. D 65 (2002) 065015 ; quant-ph 0012137 . 

[7] R. P. Feynman, Phys. Rev. 97 (1955) 660 . 

[8] K. Mano, Progr. Theor. Phys. 14 (1955) 435 . 

[9] R. Rosenfelder and A. W. Schreiber, Phys. Rev. D 53 (1996) 3337 ; 3354 ; 

A. W. Schreiber, R. Rosenfelder and C. Alexandrou, Int. J. Phys. E 5 (1996) 681 ; 
A. W. Schreiber and R. Rosenfelder, Nucl. Phys. A 601 (1996) 397 ; 
C. Alexandrou, R. Rosenfelder and A. W. Schreiber, Nucl. Phys. A 628 (1998) 427 ; 
N. Fettes and R. Rosenfelder, Few-Body Syst. 24 (1998) 1 . 



19 



[10] R. Rosenfelder and A. W. Schreiber, Eur. Phys. J. C 25 (2002) 139 . 

[11] See, e.g. M. J. Strassler, Nucl. Phys. B 385 (1992) 145 ; 

M. Reuter, M. G. Schmidt and Ch. Schubert, Ann. Phys. 259 (1997) 313 ; 
Ch. Schubert, Phys. Rep. 355 (2001), 73 . 

[12] C. Alexandrou, R. Rosenfelder and A. W. Schreiber, Phys. Rev. A 59 (1999) 1762 . 

[13] C. Alexandrou, R. Rosenfelder and A. W. Schreiber, Phys. Rev. D 62 (2000) 085009 . 

[14] R. Rosenfelder and A. W. Schreiber, to be published. 

[15] A. W. Schreiber, T. Sizer, A. G. Williams, Phys. Rev. D 58 (1998) 125014 . 

[16] D. Evans, J. W. Moffat, G. Kleppe and R. P. Woodard, Phys. Rev. D 43 (1991) 499 . 

[17] A. W. Schreiber, R. Rosenfelder and C. Alexandrou, in: Proceedings of the Workshop on Light- 
Cone QCD and Nonperturbative Hadron Physics, Adelaide (Australia) , 13- 22 December 1999, 
World Scientific (2000), p. 295 . 

[18] C. Chicone, S. M. Kopeikin, B. Mashhoon and D. G. Retzloff, Phys. Lett. A 285 (2001) 17 ; 
D. Zwilliger: Handbook of Differential Equations, Academic Press (1989), Chapter 42 . 

[19] M. W. Leen, Am. J. Phys. 62 (1994) 393 . 

[20] R. Rosenfelder and A. W. Schreiber, Phys. Lett. A 284 (2001) 63 . 

[21] V. E. Markushin, R. Rosenfelder and A. W. Schreiber, Nuovo Cim. 117 B (2002) 75 . 

[22] F. J. Dyson, Phys. Rev. 85 (1952) 631 . 

[23] A. Kizilersii, A. W. Schreiber and A. G. Williams, Phys. Lett. B 499 (2001) 261 . 

[24] D. L. Bennet, H. B. Nielsen and R. P. Woodard, Phys. Rev. D 57 (1998) 1167 ; 
R. P. Woodard, Phys. Rev. A 67 (2003) 016102 . 



20 



